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We revisit correlation effects in doped metallic zigzag carbon nanotubes by using both the one-loop 
renormalization group and non-perturbative bosonization techniques. Note that, if a nanotube is 
placed near a conducting plate, the long-range Coulomb interactions are screened and the resulting 
short-range interactions can be modelled by on-site and nearest-neighbor repulsive interactions U, 
V and V± respectively. Using both analytic and numeric means, we determine the phase diagram 
of the ground states. For U/t < 0.5 (t is the hopping strength), dynamical symmetry enlargement 
occurs and the low-energy excitations are described by the SO(6) Gross-Neveu model. However, for 
realistic material parameters U/t ~ 0(1), the charge sector decouples but there remains an enlarged 
SO(4) symmetry in the spin sector. 
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I. INTRODUCTION 

Carbon nanotubes are fascinating materials with re- 
markable mechanical and physical properties.^ Due to 
their low dimensionality, correlations and quantum fluc- 
tuations arising from electron-electron interactions can- 
not be ignored. It was beautifully demonstrated theoret- 
ically^ and later verified in experiments^! that the un- 
screened Coulomb interactions drive the armchair nan- 
otube into a Tomanaga-Luttinger liquid phase with ex- 
otic spin-charge separation. On the other hand, if the 
Coulomb interactions are screened so that only short- 
ranged interactions remain, various instabilities set irP 
and the ground state phase diagram is rather rich. 

As a carbon nanotube has only a few conducting chan- 
nels, it is not obvious how the long-ranged Coulomb inter- 
actions can be screened. One proposal for the reduction 
of long-ranged Coulomb interactions is a thick arrange- 
ment of nanotubes in either array or rope formPE Alter- 
natively, one can make use of improving nanoscale tech- 
nology and place a carbon nanotube close to a conduct- 
ing plate The induced image charges in the conduct- 
ing plate will render the Coulomb interactions U short- 
ranged such that U(x) ~ 1/x 3 at long distances. The 
possibility of realizing short-range interactions in carbon 
nanotubes motivated us to revisit their correlation ef- 
fects. 

While most theoretical investigations concentrate on 
the armchair nanotube because of its simplicity, it is often 
assumed that the low-energy physics should be the same 
for metallic zigzag carbon nanotubes, since for both hon- 
eycomb lattices the low-energy physics originates from 
the Dirac cones in the band structure. We looked into 
this issue carefully and found that the effective lattice 
model of the zigzag nanotube differs from that of the 
armchair. Ignoring the interactions momentarily, by in- 
tegrating out the gapped modest the metallic zigzag 
nanotube is described by an effective two-chain system 
(shown in Fig. [TJ . Note that the interchain hopping is ab- 
sent in this effective lattice model. As a consequence, the 




FIG. 1: (Color online) Zigzag carbon nanotube and the effec- 
tive two-chain model after integrating out the higher-energy 
modes. 



Fermi velocities remain degenerate when the nanotube is 
doped away from half filling. It is important to empha- 
size that the Fermi velocities in the armchair nanotube 
are equal to each other only at half filling and become 
different at finite doping. As we will explain later, the 
degeneracy of the Fermi velocities in the zigzag carbon 
nanotube is crucially important for dynamical symme- 
try enlargement (DSE) at finite doping. In addition to 
the subtle velocity degeneracy at finite doping, the dif- 
ferences between the zigzag and armchair nanotubes' ef- 
fective lattice models become obvious when one tries to 
write the effective interactions. It turns out that, after 
integrating out the gapped modes, the effective interac- 
tions for the zigzag are far more complicated that those 
for the armchair. 

We reanalyze the correlation effects in the effective 
two-leg lattice model with the one-loop renormaliza- 
tion group (RG) and the non-perturbative bosonization 
techniques. Due to the spreading of the wave func- 
tion, the effective interaction strength is reduced to 
U e B, V e fi — U/N x , V/N x , where N x is the number of unit 
cells around the circumference and U, V are the on-site 
and the nearest-neighbor repulsions. This reduction in 
the effective interaction strength is crucially important 
as it allows the weak-coupling RG analysis to be valid 
provided U/t < 2ttN x . 
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The one-loop RG flows determine the instabili- 
ties which may then be analyzed by non-perturbative 
bosonization (and rcfcrmionization) in order to identify 
gaps in different charge and spin sectors. From this, not 
only can the phase diagram be determined, but our nu- 
merical results also reveal dynamical symmetry enlarge- 
ment (lDSE]oflow-energy excitations above the ground 
stateP ^ 2 * 13 * 14 * 15 ! For a truly weak interaction strength 
U/t < 0.5, DSE occurs and the system is described by 
the Gross-Neveu model with an enlarged SO(6) symme- 
try. The exact solution gives eight- fold kink/antikink ex- 
citations in the low-energy limit. This differs from the 
usual behavior of two decoupled single chains which ex- 
hibit spin-charge separation but only have four-fold de- 
generate spinons. It is also not the same as the fully 
gapped two-chain systems, where spinons and holons are 
glued together by long-range gauge interactions, origi- 
nating from various gaps in the charge and spin sectors. 

We also numerically compute the anomalous scaling 
of the gap/coupling ratios associated with the enlarged 
SO(6) symmetry. For bare couplings with approximate 
SO(6) symmetry, previous analytic calculations based on 
perturbation theory predict a universal exponent of 1/2 
for the scaling functions of the gap/coupling ratios. Sur- 
prisingly, not only do we verify the predicted exponent 
numerically, we also find that the exponent is robust even 
for large deviations in the bare couplings! Finally, for re- 
alistic material parameters U/t ~ 0(1), the charge sec- 
tor decouples and ruins the SO (6) symmetry but there 
remains an enlarged SO (4) symmetry in the spin sector, 
indicating an unexpected degeneracy between the usual 
three-fold magnons and an additional neutral particle- 
hole excitation. 

The rest of the paper is organized in the following way: 
In Sec II, we derive the effective lattice model for the 
zigzag carbon nanotube at finite doping. By carefully 
taking the underlying lattice structure into account, we 
write down the effective interactions for the on-site and 
nearest-neighbor repulsions. In Sec. Ill, we show how 
the field theory description can be derived from the ef- 
fective lattice model by chiral-field decomposition. Then, 
the bare couplings in the field-theory limit can be deter- 
mined from the effective lattice model. In Sec IV, the re- 
sults of the one-loop RG are shown. From the bosoniza- 
tion technique, the phase diagram of the ground state 
is obtained. We also show how DSE emerges from the 
numerics in weak coupling. In Sec. V, we consider the 
anomalous scaling behavior rising from the finite strength 
of interactions. Furthermore, we show that the enlarged 
SO (4) symmetry in the spin sector is robust as interac- 
tions increase. In the last section, we discuss the physical 
implications and consequences of our results. 



II. EFFECTIVE LATTICE MODEL 

A carbon lattice may be represented by two regular 
triangular sublattices, offset by d = a(0, — l/v3) and 



with basis vectors a± — a(±l/2, V3/2), where a is the 
lattice constant. The Hubbard hopping Hamiltonian of 
such a carbon lattice is 

H = ^ I - t ± c\ Q (r)c 2a {r + d) 

-t c\ a (r)c 2a (r + a + +d) 

-t c\ a (r)c 2a (r + a_ +d) + /i.c.| (1) 

where the spatial summation goes over all lattice points 
r = n + a + + n_a_ with integral n±. The lattice operator 
Ci a describes the destruction of a fermion with spin a in 
the i-th sublattice. The hopping amplitude along the 
vertical (i.e., y direction) bond is t± and along the two 
zigzag bonds is t. However, in the following calculations, 
we only consider isotropic hopping with t = t± . 
The Hamiltonian for the on-site interaction is 

H u = uY,^(r)n ll (r) 1 (2) 

r . i 

where U is the on-site interaction strength and the elec- 
tron density is defined by ni a = c ia Ci a . Similarly, the 
nearest-neighbour interaction Hamiltonian is 

Hv = ^2 j Vj - n ^(r)n 2 p(r + d) 

r,a,/3 ^ 

+V n la (r)n 2l3 (r + a + + d) 

+V n la (r)n 2l3 (r + a_ (3) 

where V± is the nearest-neighbor interaction strength 
across the vertical bond and V is the nearest-neighbor 
interaction strength across the two zigzag bonds. In re- 
ality, we expect V, V± < U. 

The Hubbard model of the zigzag carbon nanotube 
may be mapped onto a two-leg ladder.^ Here the y-axis 
is defined to be in the longitudinal direction of the nan- 
otube and the x-axis is around the nanotube. Due to the 
translational invariance around the nanotube, advantage 
is taken of the fact that the momentum must be quan- 
tized in the x direction, 

fc. = ^, p = 0,±l,...,±(JV x /2). (4) 
aJM x 

The significant momenta are those which coincide with 
the Dirac points, i.e., the zeros of the energy spectrum. 

In the weak coupling limit U <C t,t±, it is natural to 
derive the effective model for the hopping first. Note that 
the energy spectrum for two-dimensional graphene in the 
tight-binding limit is 

h(k) = 2tcos(k x a/2y k « a/2 ' / * + t e -' k y a/ ^. (5) 

Therefore, the Dirac points are given by k — 
(±47r/3a, 0), (±27r/3a, ±27r/v3a). For a metallic carbon 
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nanotube, N x — 3n and the quantized momenta cut 
through the Dirac points at k x = ±27r/3a. For the nu- 
merics in this paper, we choose N x — 12 but the results 
hold for general metallic zigzag nanotubes. 

After a partial Fourier transformation in the x direc- 
tion and integrating out the gapped modes, the lattice 
operator Ci(r) is described by the fields at k x — ±2ir/3a, 



1 



E d ii(y) e 



iq{2ir /3a)x 



(6) 



q=± 



On substitution of Eq. (|6| into the hopping Hamiltonian, 

H = -t [<i(tf)<W) 
y,q=± 



+d\ 1 {y)d q2 {y l -2b) + h.c:_ 



(7) 



where the average lattice constant is b — ay3/4 with an 
offset between the two sublattices S — a^/i/12. That 
is to say, sublattice 1 consists of sites at y = 2mb and 
sublattice 2 at j/' = y + b — S. This effective hopping 
Hamiltonian describes a two leg ladder with no hopping 
across the rungs as shown in Fig. [I] Since there is no 
hopping between the two legs, the Fermi velocities of the 
two bands remain degenerate even at finite doping. The 
armchair nanotube can also be captured by a similar (but 
different) effective two-leg model, but in this case hopping 
between legs is permitted and the Fermi velocities will 
become different when doped away from half filling. 

Now let us turn to the derivation for the interacting 
Hamiltonians. From Eq. ([6]), the effective Hamiltonian 
from on-site repulsion is 



H„ = 



U 



^KiT(y) n «u(y) + n q it(y)ng U (y) 



y-Q 



U_ 



+4iT ) d 5iT (f )4u (v)d q ii (y)} 

Yi n i^(y') n in(y') + n>q2t(y')nq2i(y') 
y,q 

+4 2T {y')d m {y')d\ 2l {y')d q2i (y')} , (8) 

with n q i a = d qia d q i a and y' = y + b — 5. The notation 
q = — q is introduced for simplicity. Note that the above 
Hamiltonian contains not only the density-density inter- 
action between the two legs q = ± but also the pair-pair 
interaction between (g, q) singlets. Similarly, we can de- 
rive the effective Hamiltonian for the nearest-neighbor 
repulsion, 



2V 



EE 

y,p,q, a, (3 L 



n p i a {y)n q2 f}{y') 



■cos(2n/3)5 pq dl la (y)d qla (y)dl 20 (y')d p2f3 (y') 



N x 



EE 

V,P,Q afi 



n p i a {y)n q 2p{y' - 2b) 



+8 pq dl la {y)d qla {y)d\ 2l3 {y' - 2b)d p2 p{y' - 2b) 



(9) 



It is important to emphasize that the factor cos(27r/3) 
comes from the underlying honeycomb lattice after a par- 
tial Fourier transformation and can not be gauged away 
by shifting the second sublattice by S. Furthermore, since 
the electrons are delocalized around the nanotube, both 
the effective on-site and nearest-neightbor interactions 
are suppressed by 1/N X . 

The effective hopping Hamiltonian in Eq. ^ and the 
two interaction ones in Eqs. ([8]) and ^ now resemble an 
effective two- leg ladder. However, they are quite different 
from standard two-leg ladder Hamiltonians. The two- leg 
ladder Hamiltonian usually describes hopping and near- 
est neighbour interactions both along legs and across 
rungs, with the rungs positioned at right angles to the 
legs. However, in the zig-zag nanotube case the hopping 
only occurs along legs, and not across rungs. In addition, 
the effective interactions are far more complex than the 
ordinary density-density interactions. 

In particular, the inclusion of the next-nearest neigh- 
bor interaction really complicates the story since it re- 
flects the underlying lattice structure. In fact, if one goes 
through the same mapping, it is straightforward to show 
that the effective interactions of V and V±_ for the zigzag 
and armchair nanotubes are different. On the other hand, 
the effective interaction of U can be shown to be equiv- 
alent by a simple (but subtle) gauge transformation. 



III. FIELD-THEORY LIMIT 

To determine the ground state of the effective two-leg 
model, the combined usage of one-loop RG and bosoniza- 
tion is extremely powerful. In order to proceed, one needs 
to obtain a field theory description of the lattice model 
in the continuous limit. When considering weak inter- 
actions, a two-leg ladder is usually diagonalized so that 
the hopping Hamiltonian is written in terms of two de- 
coupled bands. However, in our case the two legs are 
already decoupled in the hopping Hamiltonian, making 
any diagonalization unnecessary. In the low-energy limit, 
the lattice operators may be linearized about the Fermi 
points and expressed in terms of chiral fields. For the 
first sublattice, the chiral- field decomposition is rather 
standard, 



d q i(y) w Vb 



ipRq(y)e ikFV + ip L q(y)e- ik ™ 



(10) 



For the second sublattice, one may expect the finite offset 
6 should give rise to some phase factor in the chiral-field 
decomposition. However, a careful analysis leads to a 
somewhat surprising result, 



d q2 {y + b-5) re Vb 



ipR q (y + b)e lk ^y+v 



+ i>Lq(y + b)e- ik ^v + V 



(11) 
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The above result states that the chiral-field decomposi- 
tion is the same as when the finite offset <5 is ignored. 
The detailed derivation can be found in Appendix |A") 

Now we are in a position to rewrite the effective two-leg 
model in terms of these chiral fields. The Hamiltonian 
density can be separated into the kinetic and interacting 
parts, H — Ho + Hi, 



Ho 
Hi 



V Yy^Rqc t id V^Rqc ~ ^l q Jd y lp Lqa ], (12) 



q , a 



C„ n iJRaa' J Laa' -cZ n ,J 



Rqq' <J Lqq' "qq'JRqq' ' J Lqq 



J Laa 1 



fqq'jRqqJhq'q' f qq > J Rqq ' J Lq' q' 



(13) 
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FIG. 2: (Color online) Phase diagram for doped zigzag car- 
bon nanotube with different V±/U and V/U. 



where v = \fZta/2 is the Fermi velocity and q,q' = 1,2 
are the chain indices in the effective model. It is worth 
mentioning again that, due to the absence of the inter- 
chain hopping in the effective two-chain system, there is 
only one Fermi velocity. The electron-electron interac- 
tions are written in terms of the SU(2) currents, 



Jpqq'iv) 

Jpqq'(y) 



l^p ga (y)<PpMy), (i4) 

^PqaiV^^lpPq'^y)- (15) 



Here P = R,L denotes the chirality of the fields and cr 
is the vector of Pauli matrices. 

Making use of operator product expansions for the cur- 
rent products, one can derive the one-loop RG equa- 
tions, dgi/dl — A^gjgk., where gi are the couplings in 
Eq. (13 1. For the effective two-chain system considered 



here, we need four couplings for the Cooper scattering in 

). 



the charge and spin sectors, c{ x ( 



) and c£ 2 ( 



and two forward scattering ones f p 2 {= $21 )■ The cou- 
pled one-loop RG equations are too complicated to allow 
analytic solutions except for some special cases. The bare 
couplings 3i(0), required for numerical integrations of RG 
flows, can be computed from the lattice interactions U, 
V and V±, 

<f xl = (U + 3V + 3V ± )/t' <? ll = (U + V- VI) A'(16) 
c p 12 = (U + W±)/t' c a l2 = (U-2V- V ± )/t'(17) 
/f a - (U + W + W ± )/t' r i2 = (U-2V- Vl)/i'(18) 

where t! — TrN x v/b — 2irN x t. Because of the 1/N X fac- 
tor, the applicability of the "weak-coupling" regime is 
extended to U/t < 2irN x . The values of t reported in the 
literaturtpE3l! range from 2.4-2.7 eV for CNTs, while 
t ~ 3 eV is typical in graphites. Although an accurate 
value of U is not yet known in nanog raphit e systems, the 
value for polyacetylene, U ~ 6-10 eVj 19 * 2 ° l might serve as 
a reasonable guess. Thus, we expect that U/t ~ 2 — 4 
in nanographite systems, which satisfies the criterion for 
"weak coupling" . 



IV. DYNAMICAL SYMMETRY 
ENLARGEMENT IN WEAK COUPLING 

To determine the phase diagram of the ground states, 
we first integrate the RG equations numerically up to the 
cutoff length scale, defined as when one of the relevant 
couplings reaches unity, gi(l c ) — 1. With N x — 12 and 
as long as U/t < 10, the couplings fall into two groups at 
I = l c : the relevant couplings reach order unity and the 
irrelevant ones remain much smaller than unity. How- 
ever, the one-loop RG alone can only determine insta- 
bilities. To truly pin down the ground state, we need to 
employ the non-perturbative bosonization and refermion- 
ization techniques. We follow the bosonization scheme 
described in Ref. [TT] The kinetic part takes the form, 
Wo = & E M PA) 2 + (9y^)% where n = p±,a± de- 
notes the total and relative charge and spin sectors re- 
spectively, with the conjugate bosonic fields # M and <p M 
describing the displacement and phase fluctuations. The 
interacting part can also be bosonized after dropping the 
irrelevant couplings in the RG flows. For instance, the 
bosonized interaction Hamiltonian of the d-w&ve super- 
conductor phase in Fig. [2] is, 



i=p+,l,2,3 

^ Bij cos 6i cos 0j , 



(19) 



i,j=l,2,3(i^j) 



where (0i,9 2 ,6 3 ) = (<Pp-> $cr+, $<>•-)• The coefficients in 
front of the gradient terms are K p +/i = /^^ih -^2/3 — 
~ c ii~^~fi2- The coefficients are symmetric with i? 12 = 
cfc, B 13 = (cf 2 + c p 12 )/2 and B 23 = -c^. 

The non-zero B^ pin the corresponding neutral bo- 
son fields in a consistent way and generate gaps in the 
$1.2,3 sectors, while the total charge mode 8 p+ remains 
gapless (described by the Luttinger liquid) and is pro- 
tected by translational invariance. With the specific 
pinned values of $1,2,3, one can show that the ground 
state exhibits quasi-long-range superconducting correla- 
tions with (i-wave symmetry while all other correlation 
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FIG. 3: (Color online) Log-log plot for anomalous scaling of 
Dij in weak coupling for U/t ranging from 10 to 1. The 
data collapses onto two universal straight lines with slope 1/2, 
as predicted by perturbation theory. 



FIG. 4: (Color online) Coupling ratios of Ki and By at the 
cutoff length scale for the N x = 12 doped zigzag carbon nan- 
otube with different bare interaction strength. The three filled 
symbols represent the relevant ratios related to the enlarged 
SO (4) symmetry in the spin sector. 



functions decay exponentially. We thus name this phase 
a d-wave superconductor (dSC). When V and V± grow 
larger, the ground state undergoes a quantum phase tran- 
sition to the charge density wave (CDW) shown in Fig. [2] 
Since the effective Hamiltonian for the CDW is rather 



similar to that in Eq. (19 1, we will not give a detailed 



derivation here. Clearly, changing the interaction profile 
allows the system to be tuned through the dSC-to-CDW 
quantum phase transition. 

The 1/N X reduction of the effective interaction 
strength gives some unexpected results. For simplicity, 
let's start with the extremely small U/t = 10~ 3 with- 
out nearest-neighbor interactions. We choose the cut- 
off length scale to be |cji(f c )| = 1- Remarkably, our 
numerics show equality among all coefficients Ki and 
Bij at the cutoff length scale. This beautiful equality 
among all coefficients signals an enlarged symmetry In 
fact, following the refcrmionization scheme developed in 
Ref. QTJ one can define six Majorana fermions i]f with 
% = 1,2, 6 to replace the three bosonic fields #1,2,3- In 
terms of these Majorana fermions, the neutral sectors of 
the Hamiltonian density in Eq. ( 19 ) becomes the exactly 



soluble Gross-Ncvcu model with the enlarged SO(6) sym- 
metry, 



Hn — 



• E 

a,b=l 



(20) 



where Q a ± describes the total/relative z-component of 
the spin between two chains and Q p - describes the 
relative charge. The six fundamental fermions cor- 
respond to the usual three-fold degenerate magnons 
and three additional particle-hole bound states and 
carry the quantum numbers (Q P -,Qa+,Qa-) = 
(±1, 0, 0), (0, ±1, 0), (0, 0, ±1). The kinks/antikinks carry 
the fractionalized quantum numbers (±|,±|,±|). The 
kinks are defined to have an even number of positive 
quantum numbers, while the antikinks have an odd num- 
ber of positive quantum numbers. Note that, since the to- 
tal charge sector decouples, all the kinks are charge neu- 
tral and can be considered to be the generalized spinons 
originally found in the single-chain system. However, 
the charge and the spin sectors are not fully decoupled 
here (as opposed to the single-chain system) because the 
kinks still carry definite Q p - — ±|. This partial spin- 
charge separation is reflected in the degeneracy of the 
spinon-likc excitations. In the single-chain system there 
are two kink/antikink excitations, so one may expect the 
degeneracy to double to four-fold in a two-chain system 
with complete spin-charge separation. However, due to 
the incomplete spin charge separation the degeneracy of 
the kinks and antikinks in our two-chain system, and 
its equivalent zigzag nanotube, is eight-fold rather than 
four- fold. 



where Gf lb — iri^rj^ are the generators of the SO(6) rota- 
tions. The excitation spectrum of the SO(6) phase from 
the exact solution consists of two distinct features: eight 
kink/antikink excitations with mass m^, and six funda- 
mental fermions with mass m\ — V2mk- Above the 
threshold E c = 2m,k is the kink/antikink continuum. 

To understand these excitations better, it is convenient 
to define the conserving charges 



(Qi - Q 2 )/2, 



)a±=Si±SZ, (21) 



V. ANOMALOUS SCALING AT FINITE 
INTERACTION STRENGTH 

The unexpected enlarged SO(6) symmetry for ex- 
tremely small U/t is beautiful but can it survive when 
the interaction strength increases? Suppose A$ denotes 
the gaps for different fundamental fermions. If the SO (6) 
symmetry is no longer exact, the gap ratios will deviate 
from unity and can be calculated from the coupling ra- 
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tios at the cutoff length scaled When DSE occurs, the 
equality among K t and B^ implies that the ratios among 
the couplings at the cutoff length scale arc universal, 
gi(l c )/9j(lc) — r i/ r ji where Ti are order-one constants 
solely depending on the enlarged symmetry. Thus, it is 
convenient to define the deviations of the coupling ratios 
as a measure for the enlarged symmetry, 



524 G^ A G^ A . 



(24) 



Di 



9ji l c) 



(22) 



where the cutoff length scale is defined by 1 0^(^)1 = 
1. If the bare couplings only deviate from the enlarged 
symmetry slightly, one can show that D^- oc (U/2nt) 1 / 2 
for the SO (6) symmetry. 15 Following a standard scaling 
argument, the gap ratios at finite interaction strength 
also show anomalous scaling, 



A* 

A; 



= i + a 



f— 

\2irN x t 



1/2 



(23) 



where are order-one constants. The universal expo- 
nent 1/2 is related to the enlarged S0(6) symmetry and 
does not depend on whether the ground state is in the 
dSC or CDW phase. 

However, these analytic results are not necessarily ap- 
plicable to our situation. For physical short-range in- 
teractions, the bare couplings deviate from the desired 
S0(6) symmetry significantly and the predictions from 
perturbation theory may become fragile. Thus, we nu- 
merically integrate the RG flows for U jt varying from 
10~ 5 to 1 and directly compute the coupling ratios at 
I = l c . The numerical results are shown in Fig. [3] It is 
truly remarkable that the scaling form with exponent 1/2 
holds true for over five orders of the interaction strength! 

Further increases in U/t > 0(1) will deviate from the 
anomalous scaling regime. Since the constants Cij are 
difficult to derive from perturbation theory, we resort 
to a numerical approach to determine at what interac- 
tion strength the SO (6) symmetry is no longer a good 
approximation. Starting from different bare interaction 
strengths, we compare the coefficients K\ and Bij at the 
cutoff length scale, as shown in Fig. [4] For U/t < 0.5, all 
coefficients deviate by less than 20% and the SO (6) sym- 
metry serves as a good approximation. However, when 
U/t > 0(1), the charge sector K\ increases rather rapidly 
and ruins the SO (6) symmetry. Note that this is also 
where the scaling form starts to break down. However, 
not all equalities are broken at the same time. Interest- 
ingly, even for U/t = 4, one can read from Fig. [1] that 
the coefficients in the spin sector remain approximately 
equal, K2 ~ K3 ~ i?23- A closer examination reveals 
that B 12 ~ B 13 also holds. 

The secret message embedded in the approximate 
equalities of the coefficients can be revealed by the 
refermionization technique. After refemionization the ef- 
fective Hamiltonian now takes the form, 



Ti-N - 92 ^ G ab g ab 

AB 



<?4 > G^hG„ 



a ,b 



7 ab^ ab 



Here the upper-case indices A, B — 1,2 while the lower- 
case ones a,b — 3,4,5,6. The SO(6) symmetry is now 
broken into the smaller SO(2)xSO(4) symmetry. The 
SO(2)~U(l) is expected because of the conservation of 
the relative charge Q p -. Even though the SO(6) sym- 
metry is destroyed in the physical regime U/t ~ 2 — 4, 
the enlarged SO (4) symmetry in the spin sector is sur- 
prising. Since the kinks/ antikinks are composed of soli- 
tons in all bosonic fields, the eight-fold degeneracy still 
holds, although the gap ratio between the fundamental 
fermions and the kinks will no longer be y2 as in the 
SO(6) regime. Meanwhile, the six-fold degeneracy of the 
fundamental fermions is broken down to four-fold. An ex- 
perimental investigation of these degeneracies in metallic 
zigzag carbon nanotubes should prove interesting. 



VI. DISCUSSIONS AND CONCLUSIONS 

In the above sections, we study the correlation ef- 
fects in doped metallic zigzag carbon nanotubes by us- 
ing both the one-loop renormalization group and the 
non-perturbative bosonization techniques. In particu- 
lar, we are interested in the ground state properties in 
the presence of short-range interactions. Note that, if 
a nanotube is placed near a conducting plate, the long- 
range Coulomb interactions are screened and the result- 
ing short-range interactions can be modelled by on-site 
and nearest-neighbor repulsive interactions U, V and V± 
respectively. 

By integrating out the gapped modes, the metallic 
zigzag nanotube can be mapped into an effective two- 
leg model. Due to the derealization of electrons around 
the tube, the effective interactions are reduced by 1/N X 
and thus can be treated as weak-coupling perturbations. 
However, since the energy of the gapped modes (which 
we integrated out) scales to 0(t/N x ), the two-leg model 
will break down when N x is too large and the multibands 
must be included. 

Using both analytic and numeric means, we determine 
the phase diagram of the ground states. For U/t < 0.5 (t 
is the hopping strength), dynamical symmetry enlarge- 
ment occurs and the low-energy excitations are described 
by the SO(6) Gross-Neveu model. However, for realistic 
material parameters U/t ~ 0(1), the charge sector de- 
couples but there remains an enlarged SO (4) symmetry 
in the spin sector. 

Finally, we thank the grant supports from the National 
Science Council in Taiwan through NSC 94-2112-M-007- 
031(HHL) and NSC 93-2112-M007-005 (HHL). 



APPENDIX A: CHIRAL FIELD EXPANSION 



In this appendix we derive the chiral field expansion of 



the fcrmion operator d g i in Eqs. ( 10 I and (111. We firstly 



7 



take a Fourier transform of the hopping Hamiltonian Eq. 
([7]) and show that the eigenvalues are E± = ±2i| cosfc&| 
for which we can define two eigenfunctions c±(k). If k 
is chosen to lie in the interval [0, 7r/6] we can write the 
original fcrmion operators in terms of the new eigenfunc- 
tions, 



1 

V2' 



d ql (k) = -^=[c q+ (k) + c q -(k)] 



d q 2{k) 



sign(fc - k F ) ikF s 

\[2 '' ' ' " ' ' '' 



[c q+ {k) - c,_(A)] (Al) 



where the Fermi momentum is kp = it /2b. 

If we take a spatial Fourier transform we need only 
retain the soft modes, which lie near the Fermi point, 



dqiiy) 



1 



V2N 

Cq+(k} 



■ (Ai 



Y^[cq+(k F +p)e 
p 

p y(k F -p)y 



i(k F +p)y 



(A2) 



where < p' < A for some cutoff A. The chiral fields are 
defined as 



■P), 
■P), 



<p < A 
-A < p < 



•p), 



< p < A 
-A < p < 0. 



(A3) 



Therefore, after taking the Fourier transform of the chiral 
fields from p to y we obtain the first part of Eq. (10 1. 
Note that y = 2mb so e lk " v = e~ lkFV . 



Similarly we may write d q 2{y + b — 5) in terms of the 
chiral fields. In this case we find that the finite offset 5 in 
the spatial coordinate is cancelled by the e lkpS which ap- 
pears in the expansion of d q 2(k) in Eq. ( Al ). Therefore, 



the lattice operator in terms of c q ±(kp ± p) is 



d q2 (y+b+s) = -L=Y,M k F+py {kF+p)iv+b) 

p 

- c q+ (k F -p)e t( - kF -P^ + ^ 



i(k F +p)(y+b) 



- Cq-(kF+p)e 

+ c q ^(k F -p)e l ^-rt (y+b *>} 



(A4) 



once all but the soft modes have been discarded. Substi- 



tuting Eq. ( A3 1 and taking the Fourier transform gives 



the second part of Eq. (Ill 
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